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We construct a new statistical physical model of poly- 
mer translocation through pore in membrane treated as the 
diffusion process across a free energy barrier. We determine 
the translocation time in terms of chain flexibility yielding 
an entropic barrier, as well as in terms of the driving mech- 
anisms such as transmembrane chemical potential difference 
and Brownian ratchets. It turns out that, while the chemical 
potential differences induce pronounced effects on transloca- 
tion due to the long-chain nature of the polymer, the ratchets 
suppress this effect and chain flexibility. 

PACS number(s): 05.40. +j, 83.10.Nn, 87.22.-q, 87.22.Fy 
(To appear in Phys. Rev. Lett.) 

The process of polymer translocation into or across 
biomembranes is a problem of considerable importance 
to a multitude of biological functions. Proteins are trans- 
ported across a cellular membrane and endoplasmic retic- 
ulum, while RNAs across a nuclear membrane after their 
synthesis There are similar macromolecular trans- 

port mechanisms in drug delivery, as well as in biotech- 
nology of gene transfer |3| where it is fundamental to un- 
derstand how DNAs can be incorporated into cells. It is 
a highly complex process with specificity involving con- 
formational changes of the translocating polymers that 
can occur in both cis and trans sides as well as inside of 
the membrane. 

Although the translocation apparatus have been sug- 
gested and examined empirically in a great variety [|4|, 
only recently there have been a few efforts to investi- 
gate quantitatively the driving force of translocation on 
physical grounds |M. Baumgartner and Skolnick || 
studied via simulation the translocation of polymer di- 
rectly through lipid bilayer, driven by the concentration 
imbalance of lipids that exists at high-curvature regions 
in membrane. On the other hand, Simon, Peskin and 
Oster(SPO) |j) considered protein translocation through 
a translocation channel or pore, and postulated that 
its driving force is random thermal motion rectified by 
'ratchets' which give rise to directional diffusion. The 
origin of this so called Brownian ratchets (BRs) is a chem- 
ical asymmetry, i.e., if specific predetermined segments of 
the protein cross the membrane, chemicals such as chap- 
erones bind on the segments to prevent their backward 
diffusion to the cis side of the membrane. 

SPO considered rigid proteins, leaving out the effects 
of three dimensional chain conformations and the asso- 
ciated flexibility and entropy. In this Letter we incorpo- 
rate these important effects by considering the flexible- 
polymer model. Along with the BR mechanism, wc 
also incorporate the more ubiquitous kind of asymme- 



try due to transmembrane chemical potential difference, 
which naturally exists in biomembranes, due to e.g., elec- 
trochemical gradients, membrane potentials and protein 
conformational changes. We aim at an analytic, quanti- 
tative theory on the basis of statistical physics of poly- 
mer and stochastic processes. To highlight these flexibil- 
ity and asymmetry effects on translocation, we consider a 
simple, but tenable model for the membrane: a rigid wall 
of negligible thickness with a pore, which is assumed to be 
small enough to allow only a single segment passage. The 
interaction between membrane and polymer segments is 
considered to be only of steric origin, i.e., the segments 
cannot cross membrane except through the pore. We 
describe the translocation dynamics as a stochastic pro- 
cess crossing the free energy barrier calculated from chain 
configuration partition function. The translocation time, 
given as the mean first passage time for this barrier cross- 
ing, is obtained from the Fokker-Planck equation that we 
formulate below. The initial targeting of nascent chain to 
the pore is regarded as a separate process and is excluded 
in this study. The controversies |^,[?J over chain confor- 
mations and chaperone functions go beyond this inves- 
tigation, which is mainly concerned with finding some 
nonspecific physical principles behind translocation. 

Free Energy Barrier of Polymer Translocation — The 
conformation of a flexible polymer during its transloca- 
tion is significantly affected by steric interaction with the 
membrane, leading to a reduction of the polymer entropy 
and increase of its free energy. We adopt, as our model, 
an ideal chain with N(^> 1) Kuhn segments each with 
length b. First consider a chain with n Kuhn segments 
with the initial segment anchored on a rigid wall intro- 
duced in yz plane. With the boundary condition(BC) 
that the other segments do not cross the surface, the 
G(r, ro; n), the probability of finding the end segment at 
r, given initial one at ro on surface, is obtained using the 
image method ||; it is given as the probability for all 
configurations in free space, the Gaussian distribution 
G (r,r ;n) = (2vrn6 2 /3)- 3 / 2 exp[-3(r - r ) 2 /(2n6 2 )], 
minus the probability for the chain crossing the surface 
G (r, -r ;n), 



G(r, r ; n) = G (r, r ; n) - G (r, -r ; n) 
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where ro = (e,0, 0), e is an arbitrarily small dis- 
tance of the anchored segment from the surface. The 
steric constraint factor of chain, given as Zs(n) — 
j x>0 G(r,r ;n)dr < 1, scales as n^ 1 / 2 . In the ab- 
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sence of the constraint the partition function is given 
by Zb(ti) ~ exp(— /3n/i), where (3 — l/kpT and \i 
is the chemical potential per segment defined by fx = 
(dF(n)/dn)T in the limit n — > oo. The F(n) is the free 
energy given from the full partition function, F(n) — 
— fc^Tlog [Zs{n)ZB(n)) — \ kpT log n+/m+const, where 
the constant term is independent of n. 



FIG. 1. Schematic figure of the configuration of a translo- 
cating polymer. 

The whole chain during translocation can be decom- 
posed into two independent end-anchored chains each in 
the opposite half spaces. For the decomposition into n 
and N — n segments as shown in Fig. I, the total free 
energy is 



Tin) = F{n) + F(N - n) 



(2) 



= -fc B Tlog [n(N ~ n)} + nAp + const, (3) 

where A/i is the excess chemical potential per segment 
of trans side relative to that of cis side. The free energy 
with AyU = has a symmetric barrier of entropic origin 
which, for a long chain, is nearly flat except near n = 1 
or n = N — I(B of Fig. 2). As also shown in Fig. 2, for 



FIG. 2. Free energy T(n) 
coordinate n. (N = 1026, A: 
/3A/x = -10/iV ) 



as a function of translocation 
/3Afi = 10/N, B: (3Au = 0, C: 



a very long chain, a minute chemical potential differ- 
ence^. g., A/i = 10 _2 fcsT) can break the barrier shape 
symmetry and its contribution can dominate the free en- 
ergy. This contribution, which does not not appear for 
a polymer in homogeneous media, can yield pronounced 
effects on a translocating polymer we shall see below. 

Stochastic Model for Translocation Dynamics — For 
the long-time scale behavior of translocation, we con- 
struct a coarse-grained description in terms of the 
translocated segment number (translocation coordinate) 
n adopted as a relevant stochastic variable and in terms 
of the associated free energy barrier. It can be treated 
as a diffusive random process, which is described by a 
Fokker-Planck equation for P(n,t), the probability dis- 
tribution of n, 
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P(n,t) = C FP (n)P(n,t), 



(4) 



where CFp{n) is the operator, CFp{n) = 
l/b 2 (d/dn)D{n) exp(-0T(n))(d/dn) exp(f3T(n)). Here, 
D(n) is the chain diffusivity during translocation. In 
the case that the D remains constant, it is given by 
D = kgT/T ~ N~ u , where T is the chain friction propor- 
tional to N". The exponent v is 1 if the hydrodynamic 
interaction between the segments is neglected (as in the 
Rouse model), and is 1/2 if it is included (as in the Zimm 
model) §. 

The mean first passage time r(n,no), which is 
defined as the time for diffusion from the coordi- 
nate no to n, is obtained by solving the equa- 
tion Q £ FP (n )T(n, n ) = -1, where C FP (n ) = 
1/b 2 exp(pT(no))(d/dno)D(n ) eiq>(-pT(n ))(d/dn ). 
To obtain the translocation time for the case that 
only the front segment in trans side is ratcheted, 
we assign the reflecting and absorbing BCs respec- 
tively at n = 1 and n — N — 1: J(n = l,t) = 
-[{D(n)/b)(d/dn + pdT/dn)P(n,t)} n=1 = 0, and 
P{n = N - l,t) = 0. Under these BCs, the translo- 
cation time, defined by r = t(N — 1, 1) is integrated to 
be 
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Let us first assume, for simplicity, that D does not 
change in the course of translocation. In case of rigid 
chain without chemical potential difference A/i, T(n) — 
const, the translocation time is simply reduced to r = 
L 2 /2D ~ L 2+v , the result for the one dimensional dif- 
fusion of a single Brownian particle. Here L = Nb is 
length of the whole chain. To incorporate the chain flex- 
ibility effect, the free energy function in Eq. || should be 
included in Eq. S, resulting in, for A/i = 0, 
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r(A/i = 0) = '—^~L 2+V 



(6) 
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While the length scaling behavior of the translocation 
time of flexible chain is same as that of rigid chain, the 
prefactor of 7r 2 /8 indicates that the chain flexibility re- 
tards translocation by 23%. This trend is opposite to 
what SPO obtained, due to the entropy effect associated 
with the three dimensional chain conformation which 
they did not include ||. The translocation time is pro- 
portional to N 2+v and, remarkably, with v = 1 this scal- 
ing behavior is identical to that of chain reptation time 
in entangled polymer systems 

If there is a nonvanishing chemical potential difference, 
the translocation time can be calculated, having the an- 
alytical expressions for limiting cases, 
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(7-a) 
(7-b) 

where /i* = TV/3A/1. When the chemical potential per 
segment is reduced on trans side, the translocation time 
as given by Eq. 7-a and Eq. 7-b encounters a crossover 
in the scaling behavior from t ~ N 2+v ~ L 2+I/ to 
t ~ N 1+v ~ L 1+l/ . As shown in Fig. 3, this crossover 
occurs around fj,* = 1 corresponding to Afi = fe^T/TV , a 
very minute chemical potential difference for a long chain. 
This remarkable sensitivity of translocation to chemical 
potential asymmetry is even enhanced for the opposite 
case of higher chemical potential on trans side. Con- 
sider, for an example, a chain with TV = 10 3 , and |A/z| = 
10~ 2 /cbT, then = 10. While, this small chemical po- 
tential difference with negative sign, speeds up the poly- 
mer translocation by the factor r(/i* = — 10)/V(/.t* = 
0) = 0.18, the one with positive sign slows it down by 
the factor of 1191. Regardless of chain flexibility, this 
extreme sensitivity, already implied in Fig. 2, is a coop- 
erative phenomenon arising from chain connectivity; the 
segments respond all hand in hand (involving the scaling 
variable (i* = N/3A[a in Eq. 7) rather than as individuals 
(involving /3A^) to a driving asymmetry. This is 



reminiscent of the cooperative effect of a slight segmental 
bias that gives rise to fast protein folding as proposed by 
Zwanzig, Szabo and Bagchi Jl2| . 

The chain diffusivity can also change during translo- 
cation. Adopting the Rouse model, Z)(n) _1 = 
N^inD^ 1 + (N - njDf 1 ) where D c and D t are the 
diffusivitics of the whole chain in cis and trans sides re- 
spectively. The effect of the change, AD = Dt — D c , on r 
can be incorporated analytically, but the result does not 
affect the dramatic effect of A[i discussed above. The 
relative insensitivity of r to AD is obvious since, while 
A/i appears exponentially, D is involved inversely in Eq. 
1 

Many-Ratchet Effect — The BR mechanism, which 
was originally suggested by SPO as a nonspecific driv- 
ing mechanism for biased diffusion, assumes fast chem- 
ical binding of chaperones on chain entering the trans 
side of membrane ||. The binding sites are assumed 
to be uniformly distributed with an interval of 8 along 
the chain. To incorporate this ratchet mechanism within 
our model for the case of instantaneous action of ratch- 
ets without dissociation, the whole space of translocation 
coordinate is divided into intervals of length S = L/M, 
where M is the number of binding sites. Then the range 
of i-th interval is [i — l)a + 1 < n < iot + 1, where 
i = 1, 2, 3, ■ • • , M, and a = N/M is the number of poly- 
mer segments in each interval. The dynamics is now con- 
secutive translocation(uni-directional diffusion) of each 
interval subject to the free energy therein, as well as to 
the BCs at both borders of the interval, reflecting BC 
at the left border and the absorbing BC at the right. 
These BCs are written as, J(n = (i — l)a + 1, t) = and 
P(n = ia. + 1, t) — 0, for all intervals. 

Assuming D is constant, the translocation time of the 
whole polymer is then r = T «i where 



Tj = 

D 



ia-\-l 



dn'e-^ n '\ (8) 



FIG. 3. Translocation time (in units of b 2 /2D , Do = ND) 
versus chain length TV for v = 1. (A: /3A/J, — 0., B: 
/3A/j = -10" 4 , C: fiAfj, = -1CT 3 , D: f3Afi = -1(T 2 , E: 
/3A/i = -10" 1 , F: pAfi = -1.0, G: 0Afi = -10. ) The 
crossover behavior from r ~ TV 3 to r ~ TV 2 occurs when 
TV is around fc s T/|A/i|. 



FIG. 4. f2(/i*, M) as a function of M for different values of 
// = Nf3Afi. ( TV = 1026, and A: // = 100., B: fi* = 10-74, 
C: fi* = 0., D: // = -10.74, E: // = -100. ) When M is 
sufficiently large, Q(fi*,M) goes to 1, the ratchet limit. 
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If the chain flexibility and A/z are neglected, i.e. J-(n) 
cons' . it is reduced to r = LS/2D — L 2 /2DM naturally, 
the reduction by the factor of 1/M compared with r(M = 
1), as given by SPO. In general the translocation time can 
be written as 



F Q(/z*,M). 



2DM"" r ~ , ~ J - (9) 
Numerical integration for translocation time gives 
Q(/z*,M) as depicted in Fig. 4, which clearly indi- 
cates that the ratchets suppress the chain flexibility, as 
well as the chemical potential difference regardless of its 
sign. Most striking is the approach of r to that of rigid 
chain(fi(/i*, M) = 1), i.e., solely the ratchet result, when 
M is very large, even with large negative value of /z*; it 
runs counter to the intuition, according to which, the neg- 
ative chemical potential difference and the ratchet mech- 
anism add up in series in speeding translocation. 

This overriding effect of many ratchets can be better 
understood by considering the Langevin equation equiv- 
alent to the Fokker-Planck equation description, 

1 dT R {n) 
b dn 

where is the Gaussian, white noise connected to T 
via fluctuation-dissipation theorem(FDT), (£(i)£(0)) = 
2TkBTS(t). Confining ourselves to the case in which the 
chain is rigid, the Tn{n) shown in Fig. 5 is the free en- 
ergy(ratchet potential) which effectively includes the BCs 
due to ratchets. For the reflecting and absorbing BCs we 
considered, the step height h is infinity. (But for general 
consideration and more realistic ratchet activities, it can 
be put to be finite. The similarity of this ratchet poten- 
tial to those employed for ratchet-driven motor proteins 
Jl3[ is remarkable.) The FDT assures the approach to 



bTh 



(10) 



FIG. 5. The free energy barrier of a rigid chain with M — 5 
ratchets and with chemical potential difference A/z. Here 
a = N/M is the number of chain segments per ratchet and 
aAzz, /i(3> a/S.fi) are the barrier heights due to asymmetries 
arising from chemical potential difference and from ratchet 
activity, respectively. 



equilibrium, that is, under the potential J-n(n) of Fig. 5, 
n undergoes the rectified diffusion to the right, regardless 
of A//. When M is very large, so that aA/z = k B Tpb* /M, 
the barrier height due to A/i, is very small, the global 
translocation dynamics become independent of the de- 
tails of the local potential barriers, yielding the t of many 
ratchets, r = L 2 /2DM, or Q(/z*,M) = 1. A calculation 
shows that this result tends to be valid for finite values 
of h larger than ksT. 

To summarize, we have investigated mechanisms af- 
fecting polymer translocation through a pore in mem- 
brane. It is found that, while chain flexibility, due to 
an entropic barrier it gives, does retard translocation, 
the ratchets speed it up, tending to reduce the flexibil- 
ity and chemical potential effects to rigid chain behav- 
ior. The transmembrane chemical potential asymmetry, 
only with minute magnitude, is found to modulate dra- 
matic changes in translocation behaviors of long poly- 
mers, which is a cooperative behavior arising from their 
chain connectivity. 
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